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Abstract 



We use a multiply- subtracted Omnes dispersion relation for the form factor in 

^ ■ B ^ n semileptonic decay, allowing the direct input of experimental and theoreti- 

^ . cal information to constrain its dependence on q^, thereby improving the precision of 

ty-^ ! the extracted value of \VuhV Apart from these inputs we use only unitarity and analyt- 

CN \ icity properties. We obtain \Vub\ = (4.02 + 0.35) x 10"^, improving the agreement with 

^ ■ the value determined from inclusive methods, and competitive in precision with them. 

^ ' 
O 

^ ' 1 Introduction 

CLiI The magnitude of the element y„/, of the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing ma- 
^ I trix plays a critical role in testing the consistency of the Standard Model of particle physics and, 
in particular, the description of CP violation. Any inconsistency could be a sign of new physics 
beyond the standard model. y„/, is currently the least well-known element of the CKM matrix and 
improvement in the precision of its determination is highly desirable and topical. 

\Vuh\ can be determined using inclusive or exclusive charmless semileptonic B decays. The in- 
clusive method has historically provided a more precise result, but recent experimental [1-4] and 
theoretical developments [5-11] are allowing the exclusive method to approach the same level of 
precision. It is important to check the compatibility or otherwise of results from the two methods, 
which currently agree only at the edge of their respective one-standard-deviation errors. 

In principle, a comparison using a calculated form factor, which contains the nonperturbative QCD 
input, at a single value of with an experimentally determined differential decay rate at the same 
would allow the extraction of \Vi,b\- In practice, experimental results are available for the dif- 
ferential decay rate integrated over q^ bins [1-4], providing shape information, while theoretical 
calculations of the form factors provide normalisation at a set of q^ values. 

Lattice QCD, originally in the quenched approximation [12-18] and more recently using dynami- 
cal simulations [8-10], provides form factor values for the high q^ region because of the limitation 
on the magnitude of spatial momentum components. Light cone sumrules (LCSR), in contrast, 
determine the form factors in the low momentum transfer region at or near = [1 1, 19-26]. 

To combine the theoretical and experimental information requires a parameterization of the rel- 
evant form factor, f^(q^), ideally based on general principles. A dispersion relation motivates 
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parameterizations by the B* pole plus a sum of effective poles (restricted and/or simplified sums 
are used in [11,27]), with a constraint imposed by the asymptotic behaviour of at large [6]. 
An alternative parameterization stems from the fact that the Bn contribution can no more than sat- 
urate the production rate of all states coupling to the uy^b current. The latter 'dispersive bound' 
was first used in this context to bound the form factors [28, 29]. More recently, it has been used to 
motivate a particular functional form which makes it easy to test consistency with the bound [5-7]. 

Here, we use a multiply-subtracted Omnes dispersion relation to obtain a parameterization of the 
form factor based only on the Mandelstam hypothesis [30] of maximum analyticity, unitarity and 
an application of Watson's theorem [31]. The latter theorem implies that has the same phase as 
the elastic nB nB scattering T-matrix in the = 1 , isospin-1/2 channel, 

= .>.,.(„„+„j^ rw = -1). (1) 

f^{s - le) T{s - le) A^'^(s) ^ ' 

The (n -1-1)- subtracted Omnes representation for f^iq^), with < s^h, reads (for more details see 
the discussion and example in the appendix of [32]): 

n n 

riq^) = m[f\sdr^'"Aoxp(ls(q';so,...,s,)]^(q^-sA (2) 

i=0 j=0 

Is(q ; so,...,Sn) = - -, (3) 

^ J.s,h (* - *o) • • • (* - s„) s-q^ 

" s s " 
ai{s) = n aiisj) = 6ij, _^Qri(^) = l. (4) 

This representation requires as input the elastic nB — > nB phase shift 6(s) plus the form factor 
values {f'^isi)} at n -l- 1 positions {sj} below the nB threshold. As the subtraction points coalesce to 
some common so, our result reduces to an expression involving the form factor and its derivatives 
at So (such a representation was used successfully to account for final state interactions in kaon 
decays [33]). The asymptotic behaviour of imposes a constraint on the subtractions (when 
more are used than needed for convergence) [34], but we keep in mind that we will apply the 
representation above only in the physical region of q^ for B ^ n decay. 

As the number of subtractions increases the integration region relevant in equation (|3]) shrinks. If 
this number is large enough, knowledge of the phase shift will be required only near threshold. 
Close to threshold, the ;?-wave phase shift behaves as 

6{s) = ni,n - p^a + • ■ • (5) 

where rih is the number of bound states in the channel (Levinson's theorem [35]), p is the nB center 
of mass momentum and a the corresponding scattering volume. In our case = 1 if we consider 
the B* as a nB bound state. Moreover, m^. is not far from Sth- We will perform a large number of 
subtractions so that approximating 6(s) x n in equation (|3]) is justified. The factor can then be 
evaluated analytically and we find an explicit formula for f^{q^) when q^ < 5th, 



(6) 



This amounts to finding an interpolating polynomial for ln[(5th - q^)f^{q^y\ passing through the 
points ln[(5th - ^- )/""(<?• )] at q]. 
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In equation Q we have assumed that has no poles. In the Omnes picture, the B* is treated as 
a bound state and is incorporated through the phase-shift integral. Since m|. is close to s^^^, the B* 
pole's influence appears in the factor 1 /(i'th - q^) in equation Going beyond the approximation 
6{s) = n, the form factor will be sensitive to the exact position of the B* pole, since the eff"ective 
range parameters (scattering volume, . . . ) will depend on mg.. 

In the following we use the explicit formula in equation ^ with four subtraction^. We have 
performed a simultaneous fit to values from unquenched lattice QCD and LCSR calculations, 
together with experimental measurements of partial branching fractions. Our main results are: 

\VJ = (4.02 + 0.35) X 10-^ Wjrm = (8.7 + l.O) x lO-\ (7) 

The 9% error for \Vub\ is competitive with the 7% error currently quoted for the determination of 
\Vub\ from inclusive semileptonic B decays. Our fitted form factor is consistent with dispersive 
constraints [5,6]. 

2 Fit Procedure 

The hadronic part of the 5° —> n'V'vi decay matrix element is parametrized by two form factors as 

(2 2\M 2 2 

in — in \ in — m 

pb + p.- 1^^] /'(^') + / " Aq') (8) 

where q'^ = {ps - p„y is the four-momentum transfer. The meson masses are nis = 5279.4 MeV 
and m„ = 139.57 MeV for 5" and tt", respectively. The physical region for the squared four- 
momentum transfer hO < < q\^^ = {mb - in„)^. If the lepton mass can be ignored (l = e or //), 
the total decay rate is given by 

with A(q^) = (ml + ml - q^f - ^m\ml. 

Results are available for partial branching fractions, over bins in q^. The tagged analyses from 
CLEO [1], Belle [3] and BaBar [4] use three bins, while BaBar's untagged analysis [2] uses five. 
CLEO and BaBar combine results for neutral and charged 5-meson decays using isospin symme- 
try, while Belle quote separate values for 5" — > n'l^vj and B* — > tjPI'^vj. For our analysis, for the 
three-bin data, we have combined the Belle charged and neutral 5-meson results and subsequently 
combined these with the CLEO and BaBar results. Since the systematic errors of the three-bin 
data are small compared to the statistical ones, we have ignored correlations in the systematic er- 
rors and combined errors in quadrature. For the five-bin BaBar data, we assumed that the quoted 
percentage systematic errors for the partial branching fractions divided by total branching fraction 
are representative for the partial branching fractions alone and, following BaBar, took them to be 
fully correlated. 

To compute partial branching fractions, we have used r^o = 1 /Fxot = (1.527 + 0.008) x 10"^^ s [37] 
for the 5*' lifetime. 

'For four subtractions, we have checked that there are negligible changes in our results if the model in [36] for the 
phase shift is used in the integral in equation (O. 
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GeV^ 




LCSR[11] 





0.258 + 0.031 


FNAL [9] 


15.87 


0.799 ± 0.058 




18.58 


1.128 + 0.086 




24.09 


3.263 + 0.324 





fr 


GeV^ 




HPQCD [8] 15.23 


0.649 + 0.063 


16.28 


0.727 + 0.064 


17.34 


0.815 + 0.065 


18.39 


0.944 + 0.066 


19.45 


1.098 + 0.067 


20.51 


1.248 + 0.097 


21.56 


1.554 + 0.156 



Table 1 Form factor inputs for \he function defined in equation dTOl ). For HPQCD and FNAL the error 
shown is statistical only: the systematic error for input value /.'" is j//.'", where i/ = 0.10 or 0.11 respectively. 
The FNAL inputs are as quoted in [5]. 

We implement the following fitting procedure. Choose a set of subtraction points spanning the 
physical range to use in the Omnes formula of equation Q. Now find the best-fit value of \Vub\ 
and the form factor at the subtraction points to match both theoretical input form factor values and 
the experimental partial branching fraction inputs. The;^f^ function for the fit is thus (this is very 
similar to thex^ minimisation used in [5]): 

11 

i,j=l 

8 

+ J] [Bf - 5,"(|y„,|, /o, /i, /2, /3)] C,\, [Bf - fo, fufi, h)] , (10) 

kj=\ 

where y^'" are input LCSR or lattice QCD values for f^{q]) and are input experimental partial 
branching fractions. Moreover, /°™"'^'*(^f,/o,/i,/2,/3) is given by equation Q with four subtrac- 
tions {ql riqj)) at (0, fo), (qiJ3,f,), (2qiJ3, /a) and (qi,,, f,). The branching fractions 5°-"- 
are calculated using /o™'^^ The fit parameters are fo, fi, fj, fi and |y„/,|, where the latter parameter 
is used when computing 50™"'='*. We have assumed that the lattice QCD form factor values have 
independent statistical uncertainties (cr,) and fully-correlated systematic errors (6,), leading to an 
11x11 covariance matrix with three diagonal blocks: the first 1x1 block is for the LCSR result 
and the subsequent blocks have the form Cjj = cr^Sij + 6,6y. The covariance matrix, Cb, for the par- 
tial branching fraction inputs is constructed similarly with three diagonal entries for the three-bin 
inputs, together with a block for the five-bin inputs. All the inputs are listed in tables [T] and [2l 

A fit to the experimental partial branching fractions alone is sufficient to determine \Vub\f'^{q^)- At 
least one input form factor value is required in order to extract a result for |y„/,|, but we have used 
a set of theoretical inputs to reduce the final error on the fitted quantities and avoid relying on a 
single theoretical calculation. 
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range 


10^4" 






GeV^ 






CLEO [1], Belle [3] 


0- 


-8 


0.410 + 0.056 


0.451 +0.041 


& BaBar [4] 


8- 


16 


0.569 + 0.065 


0.448 + 0.039 




> 


16 


0.350 + 0.058 


0.397 + 0.041 


BaBar [2] 


0- 


-5 


0.30 + 0.05 + 0.06 


0.283 + 0.030 




5- 


10 


0.32 + 0.05 + 0.03 


0.280 + 0.031 




10- 


-15 


0.23 + 0.05 + 0.03 


0.280 + 0.025 




15- 


-20 


0.27 + 0.05 + 0.02 


0.267 + 0.028 




20- 


-25 


0.26 + 0.03 + 0.04 


0.177 + 0.022 



Table 2 Experimental partial branching fraction inputs for the;^^^ function defined in equation ( fTOl ). For the 
partial branching fractions in three bins, the error shown is statistical plus systematic combined in quadra- 
ture. For the five-bin BaBar data, the statistical and systematic errors are shown. We also give branching 
fractions calculated using our fitted form factor and \Vub\- 

0.1 , , 

7 0.08 

CD 



0.06 




5 10 15 20 25 



q^l GeV^ 

Figure 1 Differential decay rate with 68% CL band (shaded) together with experimental partial branch- 
ing fractions divided by the appropriate bin-width (histograms and points). Downward triangles denote 
combined CLEO/Belle/BaBar tagged analysis results, upward triangles BaBar untagged results. 



3 Results and Discussion 



The best-fit parameters and their Gaussian correlation matrix are: 

^ 1 



rm ^ /o 
r(^Lx/3) ^ /i 

r(2qiJ3) ^ f2 

/"^(^max) = /3 



(4.02 + 0.35) X 10"^ 
0.215 ±0.024 
0.374 + 0.041 
0.938 + 0.066 
6.63 + 1.28 



-0.31 
1 



-0.86 
0.04 
1 



-0.77 
0.39 
0.67 
1 



-0.52 ^ 
-0.15 
0.65 
0.24 
1 



(11) 



The fit has;^^/dof = 1.1 for 14 degrees of freedom. 

In figure [U we show the differential decay rate calculated using our fitted form factor and \Vub\- 
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5 10 15 20 25 



Figure 2 Form factor f^{q^) with 68% CL band (shaded) together with LCSR and lattice QCD inputs 
(circles). Downward (CLEO/Belle/BaBar) and upward (BaBar) triangles show estimates for the form factors 
deduced from the experimental partial branching fractions assuming a constant over each bin and using 
our central fitted value of \VuhV 

Partial branching fractions calculated for the same bins as used experimentally are given in the 
last column of table [21 Our calculated total branching ratio turns out to be (1.3 + 0.08) x 10"'*, in 
good agreement with (1.34 + 0.08 ± 0.08) x lO"'* quoted by the Heavy Flavours Averaging Group 
(HFAG) [37]. 

In figure|2]vve show the form factor Figure [3] shows the quantity log[(5th - s)f^{q^)l Sih\ where 
the details of the fit and inputs can be better seen. Incorporating the experimental information still 
allows a fit which is perfectly consistent with the theory form factor inputs. Note that the "exper- 




5 10 15 20 25 

q^l GeV^ 

Figure 3 Same as in figure|2]but for the quantity log[(5th - s)f^{q^)l Sf^^]. 
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0.035 
0.03 



+ 0.025 
0.02 



0.015 




-0.2 -0.1 0.1 
-z{q^, to) 



0.2 



Figure 4 Same as in figure [2] but for tiie quantity P(pf'^ plotted as a function of -z.{q^, to)- Tiie dashed line 
sitting on top of the central hne is a cubic polynomial fit to Pcpf^, see text and equation (fT2l ). 



imental" points (shown by triangles) in figures [2l [3] and |4] are obtained from the partial branching 
fractions by assuming a constant form factor over the corresponding bin and are included as a 
guide for convenience. The deviation from our curves of the highest ^^-bin CLEO/Belle/BaBar 
form factor point is not significant since the form factor varies rapidly in this region and the cal- 
culated partial branching fraction agrees within errors with the experimental one (as shown in 
table [21) . 

The inclusion of experimental shape information has balanced the tendency for the LCSR point at 
= to reduce the value of \VubV To illutrate this, using only the theory inputs and comparing to 
the total branching fraction allows the fitted form factor to pass through the LCSR point and leads 
to \Vub\ = (3.73 + 0.51 + 0.16) x 10~^, where the first error comes from the fit and the second error 
is from the HFAG total branching fraction quoted above. Moreover calculated partial branching 
fractions from this fit are above experiment at low and below it at high q^. 

We have checked that our determination of /"^ is consistent with the dispersive bound. We com- 
puted P(pf^ as a function of ziq^, to), where P, (p, z and to = ^th[l - (1 - <?max/'^th)'^^] are defined 
in reference {6^. This is shown in figure HI When P(pf^ is Taylor-expanded in powers of z, the 
constraint is that the sum of squares of the expansion coefficients is bounded above by 1 . We find 
that a cubic polynomial is an excellent fit (see figure Hj) and the coefficients are, 

a,) = 0.026 + 0.002, ai = -0.037 + 0.021, a2 = -0.103 + 0.041, aa = 0.25 + 0.37. (12) 

with 2 = 0.10;!;° ^^ < 1. The errors for the a, coefficients arise from the variation of our form 
factor Monte-Carlo propagated to P(pf^ (see the bands in figure Sj). 

One may wonder how important the inclusion of the LCSR point is for the fit. Removing this 
input leads to \Vuh\ = 4.24(40) x 10'^ /+(0) = 0.166(31), so \Vub\ increases by 6%, half its error, 
while the error itself increases by 15%. Moreover, we checked that the output percentage error in 
\Vuh\ would decrease about one-eighth as fast as the percentage error on the LCSR input decreases. 

^See equations (3), (6) and the intervening text in [6]. We use nib - 4.88 GeV and niB' = 5.235 GeV. 
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Hence the LCSR input is important for its effect on the central value, but the overall error in \Vub\ 
is not much reduced. The key to the small overall error, as noted in [5-7] is to use a model- 
independent functional form with enough parameter freedom to allow the data to determine the 
form-factor shape. The Omnes form is relatively simple and is conveniently expressed in terms of 
form factor inputs at a set of values. 

We have not included possible statistical correlations within and between the HPQCD and FNAL 
lattice inputs (the lattice analysis produces statistical correlations between the form factor values at 
different q^, while both simulations are based on the same gauge field ensembles, although they use 
different heavy-quark formalisms). We modelled correlations of the statistical errors both within 
and between the HPQCD and FNAL inputs by creating a statistical error matrix 

Cstat,7 = ro-iO-j + {1- r)aj6ij 

where r is a correlation coefficient and cr, are the statistical errors on the individual inputs quoted 
by the HPQCD and FNAL groups. We added this to the block-diagonal systematic error matrix to 
create the full covariance matrix. For r = 0.25 our fit results are essentially unchanged, while for 
r = 0.81, the central value of \Vub\ moves down by one third of the original error (away from the 
inclusive determination) while the error itself grows by 10%. We conclude that these correlations 
should be included if they are known, but unless they are strong, they will not have a substantial 
effect. 

On the experimental side, we have replaced the inputs used here with partial branching fraction 
data from BaBar in 12 bins of q^ [38], for which full correlation matrices are given. We find results 
completely consistent with those given above, but do not quote them since the data in [38] are still 
preliminary. 

Applying soft coUinear effective theory (SCET) to B ^ nn decays allows a factorisation result to 
be derived which leads to a model-independent extraction of the form factor (multiplied by \Vuh\) 
atq^ = [39]. We quote the result from our fit: 

|Vjr(0) = (8.7±1.0)xlO-^ (13) 

to be compared to |Vn/,|/^(0) = (7.2 ± 1.8) x 10*"^ in [39]. In view of this, we have tried replacing 
the LCSR input a.t q^ = with the \Vub\f*{0) constraint from SCET. The resuk here, \V„h\ = 
4.24(40) X 10 \ /+(0) = 0.167(27), is completely compatible with that using the lattice inputs 
alone QV^tl = 4.24(40) x lO-\ /+(0) = 0.166(31)). The SCET and LCSR points are not really 
compatible with each other when combined separately with the lattice inputs. Not surprisingly, the 
effects are larger on /"^(O) than on \Vub\- Finally, we also tried using both LCSR and SCET inputs, 
for which the results (\Vub\ = 3.96(34) x 10"\ /+(0) = 0.210(22)) are compatible with our quoted 
values above. 

To conclude, we have presented a theoretically-based procedure to analyse exclusive B ^ n 
semileptonic decays. Starting from very general principles we propose a simple parameteriza- 
tion for the form factor equation (|6l), requiring as input only knowledge of the form factor at a 
set of points. We have used this to combine theoretical and experimental inputs, allowing a robust 
determination of |y„/,| and of the q^ dependence of the form factor itself. Our error for \Vub\ is 
reduced compared to the current exclusive world-average value, \ Vub\ = (3.80 + 0.27 + 0.47) x 10"^, 
from HFAG [37] and is competitive in precision with the inclusive world-average value, \Vuh\ = 
(4.45 ± 0.20 ± 0.26) x 10"^ [37]. Moreover we do not find a discrepancy between our exclusive 
result and the inclusive world average. 
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